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Abstract 



CN ■ Some rigorous conclusions of the Hubbard model , Kondo lattice model 

O; 

0^ , and periodic Anderson model at finite temperature are acquired employing 



the fluctuation-dissipation theorem and particle-hole transform. The main 
S ' conclusion states that for the three models, the expectation value of S^ — S^ 

will be of order N\ at any finite temperature . 
PACS number(s): 75.10.Lp,71.20.Ad,74.65.+n,74.20.Mn 
Key words: strongly correlated, ODLRO, finite temperature 



Hubbard model(HM), Kondo lattice model (KLM) and periodic Anderson model (PAM) are three 
typical strongly-correlated electrons systems under currently intensive investigations. They exhibit 



unusual thermodynamic, magnetic and transport propertities (high-Tc) ||T|-|T6|. Dispite their su- 



perficial simplicities, exact results about them are unusually difficult to obtain in more than one 
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dimensions |T^ . Fortunately, a series of rigorous results about the ground state of the three models 
have been acquired |jl|-[0. To understand the magnetic properties of the ground state of the mod- 
els, the reflection symmetry, or up-down symmetry is often imployed and this method was initially 
utilized by Lieb|l|. Having knowing some features in the ground states, one is expecting some exact 
knowledge of the model at finite temperature . Apart from the various magnetic properties of the 
models, an extremely interesting feature is that it may provide an understanding of the high- Tc su- 
perconducting supported by cuprates such as YBaCuO. To investigate this aspect, one often utilizes 
the concept of off-diagonal long-range order (ODLRO) proposed by Yang as early as thirty years 



ago[jT8l In ||T8l, Yang showed that the existence of ODLRO of fermionic systems imply Bose-Einstein 



condensation. This relationship was made more clearer recently|jT9|] |2^. It is indeed supported by the 



BCS trial state, which does not nevertherless belong to the Hilbert space of the original Hamiltonian. 
Using the symmetry of the Hubbard model and r^-pairing, Yang constructed many eigenstates of the 
Hamiltonian supporting ODLRO 0. Later, Essler et.al. showed that the ground states of a couple 
of generalized Hubbard models possess ODLRO [^ [If. Though it is generally considered that the 
Hubbard model may account for the high-Tc enigma, the Kondo lattice model and periodic Ander- 
son model are also possibly relevent because the superconducting properties are doping-dependent 
according to experiments. 

In this letter, we make use of the fluctuation- dissipation theorem to study the pseudo-spins of 
the three models at finite temperture. The prototype one-band Hubbard model on an lattice A 

Hum = J2J2 tij4aCj<^ + UJ2 ^i^^n (1) 

(ij) <^ i 

where c\^ and Cjo- are the creation and annihlation operators of the electrons with spin a =t, i at 
site i. The hopping matrix {tjj} are required to be real and symmetric. The number operators are 
riicr = cj^Cicr, while the U denotes the on-site Coulomb interaction. It is further assumed that the the 
lattice A is bipartite in the sense that it can be devided into sublattices A and B, i.e. A = A U B, 
such that tij = o whenever {ij} G A or {ij} G B. The Kondo lattice model is 

i^KLM = E E t^Ac,. + UY: n^^n^l + JY: si- ■ S? (2) 

{ij) o- i i 

where S'"'^ are the localized spins of the impurities and S'^ are the spins of the conduction electrons ( 
whose definition will be given later). This model can be clearly regarded as a doped Hubbard model. 



The periodic Anderson model is described by 

HpAM = im Ujclcj^ + UY1 ^i1^n + J2 ^f^L + '^ II(4a/ia + flcia) + UfJ2 ^Wii (3) 

(ij) o- i icr icr i 

where raf^ = fl^fia are the number opertors of the locahzed electrons. Note that we have assumed 
that the conduction electrons have also on-site Coloumb interaction in both the Kondo lattice model 
and the periodic Anderson model. When f/ = 0, the model -ffpAM is called symmetric if 2e/ + f// = 0. 
Here f/ 7^ 0, we call the model symmetric if 2e/ -\- Uf = U and we consider this case only. 

For HM and KLM the spin Sc and pseudo-spin Sc for the conduction electrons , which is equivalent 
to the r^-pairing, are defined as follows 

s: = E 4Tc.i, S; = E cliC^r, SI = \ E(cltQT - clc,,) = \{^^- N,) (4) 

iGA jgA ieA 

St = E ^{i)QiCii, S- = E 4i)4i4v SI = - E(l - ni^ - ^a) (5) 

jgA jeA ieA 

where e(i) = Iwhen i E A and —1 when i G B .Both the spin and the pseudo-spin operators 
constitute SU(2) algebra and they commute with each other, i.e. [Sc, Sc] = , so they form an 
SU(2)(g) SU(2) algebra. For HM, the total spin is S = Ei S ■ while for KLM, the total spin is 
S = Ei(S? + S^°^). It is not difficult to show that [H, S^] = [H, S^] = [H, S,] = [H, S,] = 0. Even, 
we have [H,S] = (but [H,S] ^ and [I/klm, SJ ^ ), H = Hnu, Hklm- ■ So both HM and 
KLM enjoy SU(2)(8)U(1)®U(1) symmetry. The spin is relevent to the magnetic properties while the 
pseudo-spin is relevent to superconducting. For HM, Yang and ZhangQ showed that ODLRO exists 
whenever the expectation value of S^ — S*^ is of order A^^, where N\ is the number of the sites of the 
lattice considered. 



For PAM, the spin operators and pseudo-spin operators are] 13 



s^ = E(4Tc.i + 4/.i), s- = E(ic,T + /i/.t) (6) 

i i 

s' = l E(4tc.t - 4iC^i + fhk - fhh) (7) 

s^ = E(-i)^(c.Tca - f^^h), s- = E(-i)^(4i4T - 44) (8) 

i i 

S' = lT.i^~n,-n{) (9) 



They both constitute SU(2) algebra. The main result of this letter can be states as 
Theorem For bipartite lattice A, we have for HM, KLM. and symmetric PAM that 

< S^ - S> (^, /?) =< S, > (/i, /?)cth/5/i(^ - /i) (10) 

where < O > (/x,/5) = TrOexp(-/3K)/Tr exp(-/?K), K = H - ^iN,N = Y^i^nta for HM and KLM 
or J^iai^ia + n{„) for PAM. 



Proof. Consider the double-time Green function ^ S | S+ ^u;. The evolution equation of it is|T4 

c^ < S- I S+ >^=< [S-, S+] > + < [S-, K] I S+ >^ (11) 

It can be calculated directly that 

[S-,i^] = (2^-f/)S- (12) 

Accordingly, we have 

« S- I S+ »^= ^—— < S^ > (13) 

uj + U — Zjj, 

Therefore, the expectation value < S"'"S~ > can be obtained by virtue of the fluctuation-dissipation 
theorem 

< S+S- >= ^ r « §- I S" »"^- - « §' I §" »"-'d. (14) 

27r j-oo exp phoj — 1 

Using lim^^o ^ = ^^ =F «vr(5(x), we have 

< S+S~ >= -2 < S^ > — -^ ^ (15) 

exp/?;i(2/i-t/)-l ^ ' 

On the other hand, we have S^ — S^ = |(S+S~ + S~S+) and [S+, S"] = 2S2. Hence 

< S^ - S^ >=< S+S- > - < S, > (16) 

Using eq.(9), one can readily reach the conclusion. Q.E.D. 

Eq(15) can be employed to draw some conclusion on the electron desities in the three models. 
Since 5^ = ^kA^a - A^) where for HM and KLM, k = 1 while for PAM, k = 2 ,we have 

where pe '■= N/N\. From the finiteness of the l.h.s. of eq(17) we have the following generalization of 

the lemma in p2| 

Corollary 1 Under the same assumptions of the theorem, for all the three models, we have:pe > k 



if 2/i > U; Pe = n if 2fi = U and pe < n if 2p < U. 

For PAM, we can write S = Sc + S/. Since for a finite PAM system, we alway have 

< Se ■ S/ >=< S'J} >= (18) 

because the bra and ket do not have equal number of creation and annihilation operators of the same 
kindred. So we have the following 
Corolary 2 For PAM, 

<Sl- Sf >=< S, > cthph{^ - p)~ < S) - Sf >> (19) 

2 J J 

Since fi is the chemical potential (for PAM the conduction electrons and the localized electrons 
must share a common chemical potential because it is the total electron number instead each of the 
respect species is conserved), i.e. the fermi energy in the free fermion case, which is an intensive 
quantity, it must be a function of the particle density Pe = ^ and the temperature /3 and the model 
parameters, t, U, ...,/i = p{pe,P,t, U, ..., ) where A^e is the number of electrons accommodated in the 
lattice. Therefore, at any finite temperature,/?, which is not a root of y — p{pe, P,t,U, ...,) = 
the expectation value < S^ — S^ > can not be of order N^ as Aa -^ oo while keeping pe constant 
according to theorem 1. What about at the roots of ^ — p{pe,P,t, U, ..., ) = 0? For the HM on a 
bipartite lattice, using the complete particle-hole transform [P5[] 

PQtP^^ = e(04, PC^P-' = e{i)cl (20) 

For KLM, the complete particle-hole transform is defined by (19) together with 

PS±iocp-i ^ ^sf""", PS["^P-^ = -Si"" (21) 

i.e., the particle-hole transform effects a rotation of the spins . As to for PAM, the transform is 
defined by (19) together with 

Pf,^p-' = -eii)fl, Pflp-' = -e(z)/,., (22) 

One has for all the three models 

PKP-^ = H + k{U- 2p)Na + {p- U)N (23) 

Then since 



Therefore we always have < A^ >= kN^ whenever fi = U/2 no matter at what temperature. So 
fi = U/2 can not determine the temperature. (This is why we often say the system is at half-filhng 
when fi = U/2. for HM). We have accordingly < Sz >= in this case and the r.h.s. of (10) is in fact 
0/0, which should be determined by the limit lim^^(//2. Since 

1 d 



we have 



^^ < N >=< N^ >- <N>^ (25) 

p Ofl 

lim < S^ - S^ >= -{< N^ > - < N >2) (26) 

fj.^U/2 2 

Since for an ideal gas, the r.h.s. of (25), which is just the fluctuation squared of particle number, is 
of order 0{N), it is quite reasonable to assume the in the three models here, it is also of the order 
0{N). Hence we have therefore a stronger conclusion than that in|]2^. 



Corollary 3. Under the same assumptions of the theorem, the l.h.s. of (10) can he at most of order 
0{N\) at any finite temperature for HM, KLM and PAM. 

As in ref.||2^, eq(23) can be used to obtain a symmetry of the l.h.s. of eq(lO). Since eq(24) states 
that 

Pe(/i,/5)=2«:-pe(f/-/i,/3) (27) 

we have immediately from (17) that 

^ < S+S~ > (p, /3) = i- < ~S^~S- >iU-fx, (3) exp{/3^(t/ - 2/i)} (28) 

There exists lot of forms of the l.h.s. of (17) satisfying this symmetry. For instance, C(/3, U)e~^^^^~^^'^^ 
and C(/5, U)Y—:^jjjr:2ii) (where C{P,fi) denotes come function) . So the functional form can not be 
determined uniquely without further constraints. 

By virtue of particle-hole transform (for PAM, e is to be repalce by (—1)*) 

TQtT-i=e(z)4, Tc,^T-' = Ci^ (29) 

T/.tT-i = -e(*)4, Tfi{r-^ = f,^ (30) 

under which the spin and pseudo-spin operators transform into each other |0. 

T(5+, 5", Sz)T-^ = (5+, 5-, 5,), T(5+, ~S~, Sz)T-' = (5+, S', S,) (31) 

One can obtain the knowledge of spin of the transformed model from the known knowledge of pseudo- 
spin of a given modeip3. From the theorem we know that at half filling, < Sz >= 0, therefore we 



always have < S^ — 5^ >= at any P ^ Pc, this agrees with the Corolary 2 in [^. It was shown 
that [^ I at half-filhng for HM, no ODLRO is exhibited for on-site electron pairs in the translational 



invariant case for both positive and negative U. Yet, away from half- filling, it seems that theorem 
1 disagrees with the theorem 1 of since that theorem asserts that for negative U Hubbard model 
and some special pe, the expectation value < S^ — 5*^ > at ground state can be of order N^. But the 
special pe was given ad hoc and was not determined by the grand canonical equilibrium. 

As a by-product of the complete particle-hole transform, we have also the following 
Corollary 4. For the three models considered under the same assumptions as in the theorem, at 
half-filling, i.e. /i = U/2, one has < S^^z >= 
Proo/ Since under the complete particle-hole transform, PSx,zP^^ = —Sx,z, we have 

< Sx,z >= ^TTSx,ze-^'' = ^Tt{-Sx,z) exp{-(3[H + k{U - 2/x)iVA + (/i - U)N]} (32) 

Therefore at /i = U/2, we have 

< Sx,z >=-< Sx,z > (33) 

Hence one can obtain the assertion. 
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